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Abstract. Solid 4 He is viewed as a nearly perfect Debye solid. Yet, recent calorimetry 
measurements by the PSU group (J. Low Temp. Phys. 138, 853 (2005) and Nature 449, 
1025 (2007)) indicate that at low temperatures the specific heat has both cubic and linear 
contributions. These features appear in the same temperature range where measurements of 
the torsional oscillator period suggest a supersolid transition. We analyze the specific heat 
and compare the measured with the estimated entropy for a proposed supersolid transition 
with 1% superfluid fraction and find that the observed entropy is too small. We suggest that 
the low-temperature linear term in the specific heat is due to a glassy state that develops 
at low temperatures and is caused by a distribution of tunneling systems in the crystal. We 
propose that dislocation related defects produce those tunneling systems. Further, we argue 
that the reported putative mass decoupling, that means an increase in the oscillator frequency, 
is consistent with a glass-like transition. The glass scenario offers an alternative interpretation 
of the torsional oscillator experiments in contrast to the supersolid scenario of nonclassical 
rotational inertia. 



1. Introduction 

The discovery of an anomalous signature in torsional oscillator measurements of solid 4 He has 
re- ignited the search for the enigmatic quantum state of supersolidity [lj . Several groups using 
torsional oscillators report a drop in the resonance period around 0.15 K, with details depending 
on the amount of 3 He impurity concentration present in 4 He [U El El IH E] • The observed signal 
shows hysteresis effects depending on cooling and warming procedures or rim velocities, as well 
as annealing dependence. At this point it is far from clear if the observed phenomenon is due 
to the onset of supersolidity or the interplay of a dislocation glass or both. Very recent Monte 
Carlo simulations suggest that a quasi-one-dimensional supersolid can occur along the axis of 
screw dislocations [6J. However, it is difficult to imagine how such a network of one-dimensional 
supersolid phases is the origin for reports of as much as ~ 20% of supersolid fraction as seen in 
recent torsional oscillator experiments [?]. In this work, irrespective of the supersolid or glassy 
origin of the observed anomalies, we assume that only a small fraction of the perfect crystal 
undergoes a phase transition. 

So far x-ray and neutron scattering experiments have failed to observe any evidence of a 
supersolid condensate fraction or changes in the lattice parameter, Debye-Waller factor, or 



phonon dispersions between 55 mK and 500 mK, where the transition is reported to occur 
O [U EJ [ini E] • On the other side, specific heat [12] and elastic shear modulus [13] measurements 
exhibit anomalies in the same temperature range as torsional oscillators. 

Here we discuss how a small fraction of the sample, undergoing a glassy transition around 0.12 
K, can account for the observed linear T dependence in the specific heat of an otherwise perfect 
Debye solid. We emphasize that the physical picture invoked is not that of a conventional glass, 
where the entire sample undergoes a glass transition. Rather a small subsystem of the crystal, 
e.g., dislocations, exhibits glassy characteristics at low temperatures [14J, similar to two-level 
systems. Our glassy interpretation is consistent with the reported drop in the oscillator period 
and concomitant peak in the dissipation. The microscopic nature of the proposed glassy state 
in solid 4 He is still elusive and outside the scope of our phenomenological theories for torsional 
oscillator and specific heat data. Although lacking a microscopic model, we can outline some 
general properties attributable to this state. We propose that there is a structural component 
of the crystal (e.g., dislocations, grain boundaries or dislocation cores, etc.) that has glass-like 
features and exhibits a freeze-out of dynamics at lowest temperatures. This glassy component 
comprises a small fraction of the sample with characteristic time dynamics. In addition strong 
coupling between a glassy and possible supersolid component might be present in these crystals. 
For that one would need more detailed experiments exploring the crystal structure in order to 
address the physical nature of the glassy state at lowest temperatures. Details of the glass model 
and data analysis can be found in Refs. \15\ I16j . where it was assumed that the backaction of 
the glassy subsystem onto the torsional oscillator is small. 



2. Entropy Dilemma 

The entropy is a state function of a thermodynamic system in equilibrium and records the 
number of excited states. Thus entropy and specific heat are the preferred quantities for 
characterizing bulk phase transitions. Assuming a noninteracting Bose-Einstein condensate 
(BEC) in three dimensions with a parabolic band, the specific heat is Cbec{T) = 

^f/jf? R (T/T c ) 3 / 2 for T < T c . Correspondingly, the entropy at T c is universally Sbec(T c ) = 

Jq c (ITCbec/T = ^Mr(T/T c )V 2 ~ 5/4 fl, with gas constant R = 8.314 J/(mol K) and 

Riemann's £ function. For comparison, even the entropy of strongly correlated superfluid 4 He 
at the A point (T\ = 1.8 K and P\ = 26 bar) is of the same order, namely S\ « 0.55 R. On the 
other side, a glass comprised of a distribution of two-level systems (TLS) has a low temperature 
specific heat Ctls{T) — vR(T/T c ) for T < T c , where v is proportional to the fraction of the 
TLS in the sample. Since only a small fraction of the sample undergoes a transition into a BEC 
or glass phase, we expect that the measured low temperature specific heat for T < T c is either 
of the form C = A BEC T 3 / 2 + B Debye T 3 or C = A glass T + B Debye T 3 [TS]. 

Figures 1 through 3 show the specific heat recently measured by Lin et al. [12] on ultra- 
pure solid 4 He contaminated with 1 ppb of 3 He, which is qualitatively consistent with earlier 
measurements at higher 3 He concentrations by the PSU group [12J. Indeed, Fig. 1 shows a 
clear deviation from perfect Debye behavior below 0.12 K, with the linear-T term visible as 
a finite intercept in Fig. 2. Note, that for a BEC the C/T — > with T — > 0. Fig. 3 shows 
the difference AC/T = (C — Cr> e by e )/T, which is a measure of the entropy. We estimate that 
S(T C ) ~ 3.3 • 10~ 6 R. This value is nearly four orders of magnitude smaller than is expected 
for 1% of solid 4 He transforming into a supersolid. Therefore, neither the low-T behavior of 
C nor the excess entropy associated with a phase transition are consistent with a BEC of a 
noninteracting or strongly interacting system of vacancies. However, at this point we cannot 
rule out that below T c a dominant glass and subdominant BEC state coexist. 
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Figure 1. C/T 3 of ultra-pure (1 ppb of 
3 He) solid 4 He from Lin et al. [12]. The non- 
Debye behavior below T c ~ 0.12 K signals the 
freezing out of two-level systems. 



Figure 2. C/T of ultra-pure solid 4 He. 
Below T c , C/T deviates from the behavior of 
a perfect Debye solid with a nonzero intercept 
indicative of a small glassy component [y ~ 
3.3 • 10" 6 ). 
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Figure 3. AC/T of ultra-pure solid 4 He. 
The estimated entropy associated with the 
transition is S(T C ) 3.3 • 10~ 6 R. This 
entropy change is six orders of magnitude 
less than that for an ideal BEC transition 
Sbec(T c ) ~ 1-25 R. Thus it cannot explain 
the putative supersolid transition, even 
when assuming that only 1% of solid 4 He 
transform. 



3. Torsional Oscillator Response 

A torsional oscillator (TO) measures the susceptibility of the angular response, it does not 
directly monitor the moment of inertia of the supersolid. Hence, a change in the response 
function (period or dissipation) vs temperature is not sufficient to deduce a change in the moment 
of inertia. In linear response the angular susceptibility of an underdamped externally driven TO 
in the presence of a small backaction component g is = Xq 1 ~~ 9 [16J. In Fourier space the 
bare susceptibility of the empty cell is 

Xq T) = a- i~f osc uj - I OS cU 2 . (1) 

Here a is the restoring coefficient, 7 osc is the dissipative coefficient and I osc is the moment of 
inertia of the TO apparatus. The function g{uj; T) captures the response of added solid 4 He 
with a supersolid or glassy component that acts like a generalized polarization in a dielectric 
medium. For a solid with a supersolid subsystem one expects the backaction to be 

g ss {u; T) = i-fHe^ + I ss (T)uj 2 , (2) 

while for a solid with a glassy subsystem it is 

g g l(u) T) = i-/HeU + Igiu 2 + 5o(l - iw)" 1 . (3) 



The significant difference between both scenarios is that for a supersolid the moment of inertia 
I SS (T) is T dependent and accounts for the change in the oscillator's period, while for a glass 
I g i is T independent (to leading order) and the T behavior originates from the relaxation time 
t(T) = To exp(A/ksT) in g g i(u;T). In this expression and in our fits we assumed a simplified 
Vogel-Fulcher-Tamman relaxation time for glasses with an activated behavior, as discussed in 
more detail in Ref. |16j . 

In figures 4 and 5 we show the difference between the change in period and dissipation due 
to a glass transition vs a supersolid transition in parts of the crystal. Comparing our model 
calculations with TO measurements [2] we can reproduce both period and dissipation changes if 
we assume a glassy transition, while this is not possible for a simple supersolid transition that is 
due to a change in i# e (T) alone. The glass fit is discussed in detail in Ref. [16J. For the limiting 
case of a single relaxation mode, that means a glass exponent of (3 = 1, we were able to derive 
the simplified fit functions for period P, 

P-\T) = /o - Bf^[l + (2nf r(T)) 2 r\ (4) 

and dissipation (inverse of the quality factor Q), 

Q-\T) = Ar(T)[l + (2^/ r(T)) 2 ]- 1 + Q~\ (5) 

with frequency /o = 1/P(0). In Ref. [16] we assumed that fit parameters A and B are 
independent, however, they obey the constraint A/B ~ 47t//q I18| . This leads to 
slightly degraded fits of both period and dissipation compared to the unconstrained fits in 
Ref. [16j, though still in fair agreement with experiment. A key consequence for a glass is 
that the dissipation change, AQ^ 1 = Q^ 1 — , is given by the frequency change, AP _1 = 
— P(0) -1 |, according to AQ~ l ~ Attt AP _1 , that means, that the maximum change is 
AQ~ l ~ 2AP/P. This result is equivalent to Ref. |18j and similar to a phenomenological model 
by Huse and Khandker, up to a factor of two, who obtained AQ^ 1 ~ AP/P [19J. 

On the other hand, for the supersolid fit, we assumed the following functional dependence 
for the moment of inertia to describe the measured drop in period (note that the precise shape 
and form is irrelevant for the ensuing discussion), 

I(T) = J (l-0.5e(l- tanh[(T- 0.165 K)/0.04K])). (6) 

In order to fit the period change measured by Rittner and Reppy [2j we assumed that a small 
fraction e ~ 3.6 • 10~ 5 of the entire oscillator moment of inertia, I = I osc + Iue-, transforms 
into a supersolid. This is consistent with 0.1 — 1% of nonclassical rotational inertia fraction 
(NCRIF). Fig. 5 demonstrates that the putative supersolid transition seen in Fig. 4 does not 
account for the simultaneously reported peak in the oscillator's dissipation. It would require 
a highly unrealistic temperature dependence of the dissipative coefficient j(T) oc Q~ l (T). We 
used equations (11) and (16) of Ref. [16J to calculate the period P(T) and quality factor Q(T) 
of the supersolid with 

Q SS (T) = ^^sJl + jr-(P c /P{T))\ (7) 

and P c = (2-7T7/a) 2 < P(T). It follows that to leading order the maximum dissipation of a 
uniform supersolid is given by AQ SS /Q SS ~ AP/P or equivalently AQj*/Qj* ~ AP/P ~ 10~ 5 . 
Thus for a supersolid one expects only a very small change in dissipation Qjg, as seen in Fig. 5. 
This is in clear contrast to the glassy scenario where the maximum change is of order unity, 
AQ/Q ~ 2Q OQ AP/P ~ 1. 
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Figure 4. Change in resonance period P [2]. 
Both glass and supersolid model can explain 
the observed change. Fit parameters for glass 
are A = 0.0390s" 1 , f = 184.23119 Hz, 
= 11.33 • 10" 6 , t = 0.360 //s, A = 
1.135 K. 
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Figure 5. Change in dissipation (inverse 
quality factor Q) \2\. Only the glass model 
can explain the pronounced dissipation peak. 
Note that for supersolid scenario the small 
change in dissipation is not visible on this 
scale. 



4. Conclusions 

Our analysis of both the specific heat and the TO measurements for putative supersolid 4 He 
reveals the existence of a dilemma between the amount of the excess entropy and dynamic 
response of TO experiments. We argue that this dilemma can be overcome by invoking a glass- 
like transition below ~ 0.12 K in ultra-pure solid 4 He. Our phenomenological treatment of 
the specific heat and TO response does not allow us to identify the microscopic nature of this 
anomaly, but we speculate that parts of the solid freeze out and undergo a glass rather than 
a BEC transition. This glassy component may coexist with a possible supersolid component, 
which may or may not have been present in the samples analyzed. Thus it is reasonable to 
assume that dislocations give rise to two-level systems and form a glassy state, which dominantly 
contribute to the observed features. At this point, more dynamic studies are needed to resolve 
the physical nature of the glass-like phase and if a glass and supersolid state can coexist in the 
low temperature phase. 
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